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Abstract 

A general theorem on the GBDT version of the Backlund-Darboux 
transformation for systems rationally depending on the spectral pa- 
rameter is treated and its applications to nonlinear equations are 
given. Explicit solutions of direct and inverse problems for Dirac- 
type systems, including systems with singularities, and for the system 
auxiliary to the N-wave equation are reviewed. New results on ex- 
plicit construction of the wave functions for radial Dirac equation are 
obtained. 
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1 Introduction 

The Backlund-Darboux transformation (BDT) is well-known in the spectral 
theory and integrable nonlinear equations (see, for instance, [3,7,9,13-15,19, 
25,33,41-43,46,48,50,84] and references therein). BDT transforms initial 
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equation or system into another one from the same class and transforms also 
solutions of the initial equation into solutions of the transformed one. 

In this paper we review results on the GBDT version of the Backlund- 
Darboux transformation. GBDT works for the matrix and scalar cases, gives 
explicit expression for the iterated Darboux matrix in terms of a transfer 
matrix function, minimizes the order of the matrix that has to be inverted 
in the BDT approach. Thus, GBDT is a convenient tool to construct wave 
functions and explicit solutions of the nonlinear wave equations as well as 
to solve various direct and inverse problems. GBDT and its applications 
were treated or included as important examples in the papers [22, 23, 47, 
54-58,60,62,63,65-67,69-72,74] (sec also [27-32,36]). Here we consider 
self-adjoint and skew-self- adjoint Dirac-typc systems including the singular 
case corresponding to soliton-positon interaction and solve direct and inverse 
problems. We solve also direct and inverse problems for the system auxiliary 
to the A^-wave equation, construct explicit solutions of the iV-wave equation 
and obtain evolution of the corresponding Weyl function. The results on 
the radial Dirac equation are new and we treat them in greater detail. We 
consider also in a detailed way a general Theorem 3.1 on the GBDT for 
systems rationally depending on the spectral parameter and its applications. 
(This result was earlier published in [56].) 

First, let us illustrate BDT by the oldest and most popular example, that 
is, by the Sturm- Liouville equation 



where v{x) = v{x), A = A. Assume that z{x) = l^x) satisfies (1.1), when 
A = c, that is, —Zxx + vz = cz {zxx '■= ^^)- Then one can rewrite (1.1) in 
the form 



dx'^ 



y{x, A) + v{x)y{x, A) = Xy{x, A), 



(1.1) 





where A and A* are first order differential expressions: 




Transformed equation is given by the formula 




(1.3) 
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It easy to see that (1.3) is again Sturm- Liouville equation, but potential v is 
transformed into v — v — 2 (^—^ ■ Notice further that (^AA* + (c — A)/^ A — 

a(^A*A+ (c — A)/^ Hence, it follows that if y{x,X) satisfies (1.2), then 
y Ay satisfies (1.3). Fundamental solutions of the transformed equa- 
tions can be constructed in this way. Under rather weak conditions the 
spectra of operators A* A and A A* may differ only at zero, and so under 
certain conditions the spectra of Sturm-Liouville operators L and L associ- 
ated with differential expressions — ^ + v and — ^ + v may differ only at 
c. The Backlund-Darboux (and related commutation) methods of inserting 
and removing eigenvalues of Sturm-Liouville operators historically go back to 
Backlund, Darboux, and Jacobi [7,19,34] with decisive later contributions by 
Crum, Deift, and Gesztesy [18,20,24]. (See also [26] and a detailed account 
in Appendix G in [25].) 

One can apply BDT again to the already transformed equation (1.3) and 
so on (iterated BDT) . There is also somewhat more complicated binary BDT 
(see [1,46]). It proves that, if v satisfies nonlinear intcgrablc equation, then 
V often satisfies it too, and so BDT is used to construct solutions of the 
nonlinear equations. BDT proves especially useful for the construction of 
the explicit solutions, starting from the trivial initial system (for instance, 
t; = in (1.1)). 

Elementary Backlund-Darboux transformations for Dirac type and more 
general AKNS systems one can find, for instance, in [13,40]. Given first 
order initial and transformed systems Ux — G{x, X)u and Ux — G{x, X)u, 
their solutions are connected via so called Darboux matrix w such that Wx = 
Gw — wG. Here G, G and w are m x m matrix functions (m > 0). Clearly, 
if u satisfies some initial system Ux = Gu, then wu satisfies the transformed 
one Ux = Gu. Darboux matrix or gauge transformation is of great interest 
in this theory (see [15,40,42,46,51,78,80] and references therein). 

We assume that the fundamental solutions u of the considered systems 
are always normalized by the condition 

u{0,X)=Im, (1.4) 

where I„i is the mxm identity matrix. Then the fundamental solution of the 
transformed system is given by the formula u{x, A) = w{x, X)u{x, A)'u;(0, A)~^, 
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where u is the fundamental solution of the initial system. 

We shall consider the GBDT version of the Backlund-Darboux transfor- 
mation, where the Darboux matrix is represented in the form of the transfer 
matrix function w{x, A) = wa{x, A). Transfer matrix function corresponding 
to the S'-node (transfer matrix function in the L. Sakhnovich form) is given 
by the equality 

WAix) = I -n;s-\Ai- xiy^Ui, (1.5) 

where the matrices Ai, A2, S, Hi and 112 form an 5" — node, that is, satisfy 
the matrix (operator) identity 

AiS - SA2 = UiU*. (1.6) 

Here S and Ak {k = 1, 2) are n x n matrices and Hk {k = 1, 2) are n x m 
matrices for some integer n > 0. Sometimes we use notations Ai — A and 
Hi = n. The function wa is the transfer matrix function for system 

-^z ^ Aiz + Uiu, y^'n.lS''^z + u 
dx 

(see Introductions in [75,76]). At the same time wa is a generalization of the 
Livsic-Brodskii characteristic matrix function. 

Matrix functions Ilk{x) are introduced in terms of the generalized eigen- 
function (with the generalized matrix eigenvalue) of the initial and dual to 
initial systems (see Remark 2.4 and formulas (3.3) and (3.4) for details). 
Recall that eigenfunctions are essential for the classical Backlund-Darboux 
transformation. Matrix function S{x) is expressed directly via Ilk{x), and 
in many cases one can use for this purpose the identity (1.6). Notice that 
another type of matrix identities have been successfully used for the con- 
struction of the explicit solutions of nonlinear equations in [44]. Further 
developments of the Marchenko scheme are given in [11,12,77]. For applica- 
tion of the matrix identities to the construction of solitons see also [38]. 

Our approach grants additionally explicit expression for the Darboux 
matrix, allows to avoid the stages of the construction of the high order matrix 
solutions of the nonlinear equation and their reduction to the required order, 
and minimizes the order of the matrix that has to be inverted (matrix S in 
our case). 



4 



GBDT for the self-adjoint and skew-self-adjoint Dirac-type systems, for 
the system auxiliary to the iV-wave equation, GBDT for the A^-wave equa- 
tion itself and for nonlinear Schrodinger equation are studied in the next 
Section "Preliminaries". A much more general Theorem 3.1 on GBDT and 
some applications are given in Section 3. Radial Dirac equation is treated in 
Section 4. Various direct and inverse problems are solved in Section 5. 

By X and we denote intervals on the real axis, by diag{(ii, ^2, • • •} 
we denote diagonal matrix with the entries on the diagonal, by ^] we 
denote the commutator — ^D, and col means column. By the neighbour- 
hood of zero we mean the neighbourhood of the form (0, e) or [0, e). As 
usual we denote by Z the integers, by R the real axis, by C the complex 
plain, by R_|_ the positive semi-axis, and by C+ (C+) the open (closed) up- 
per semi-plane. By arg(a) we denote the argument of a G C. We always 
assume that X]j=A; c^r = when k > r. The notation r f means that r is a 
nondecreasing function. The spectrum of an operator A is denoted by cr(A). 
Matrices j and J have the form 



J = 



' Ip " 


, J = 


' 


Ip' 


. -Ip _ 




_ 



(1.7) 



2 Preliminaries 

In this section we apply GBDT to the self-adjoint and skew-self-adjoint Dirac- 
type systems and to the system auxiliary to the N-wave equation. We treat 
these systems on some interval I that contains 0. Applications to the A'^- 
wave and nonlinear Schrodinger equations are obtained. In this way we can 
show how GBDT works before the formulation of more general results. 



2.1 Gauge transformation of the Dirac-type system 

The self-adjoint Dirac-type system has the form 

^u{x,X)^i{Xj+jV{x))u{x,X), (2.1) 
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where u is a,n m x m matrix function, m — 2p, 



-/„ 



V 



V 
V* 



(2.2) 



and the potential v is a p x p matrix function. The skew-self-adjoint Dirac- 
type system has the form 



d 

dx 



u 



(x, A) = (iXj + iV{x))u{x,\). 



(2.3) 



Dirac-type systems are also called Dirac, Zakharov-Shabat or AKNS systems. 
System (2.1) can be rewritten as: 



— M(a;, A) + {Xqi + qo{x))u{x, A) = 0, 



Qi = -ij, Qoix) = -ijV{x). (2.4) 



GBDT for system (2.4) is generated by an integer n > and three matrices, 
that is, by n X n matrices A and S{0) — S{0)* and by n x m matrix n(0). 
It is required that these matrices satisfy the matrix identity 



AS{o) - s{o)A* = m{o)jU{oy 



(2.5) 



Fix n and parameter matrices A, S{0), and n(0). Then we define an n x m 
matrix function Il{x) by its value n(0) at x — and by the linear differential 
equation 

n,(a;) = AU{x)qi + Il{x)qoix). (2.6) 

Matrix function S{x) is easily recovered from its value <S'(0) and from the 
expression for its derivative: 

= nn*. (2.7) 

Moreover, as ^(0) = S{0)* and = S* we have S{x) = S{x)*. 

Recall that we consider systems and functions on the interval X{i.e., 
X E T), and that we assume G X. We choose point 0, where n(a;) and S{x) 
are fixed, for convenience. (It could be any other point from I.) In (2.7) and 
in some of the following formulas we omit the variables in the notations. 

Prom (2.6) and (2.7) we get 



{AS - SA' 



Aun* - IlIl*A* 
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{inju% = (mju*) - {mjuy = Ann* - nn*^*, i.e., 

{AS-SA%^{mjU%. (2.8) 
Formulas (2.5) and (2.8) imply identity 

AS{x) - S{x)A* = m{x)jll{xy. (2.9) 

That is, matrices A, S{x), and n(x) satisfy identity (1.6), where 

Ai = Al = A, ni = n, Ul = ijU*, (2.10) 

and we say that A, S{x), and Il{x) form an 5'-node. Using (2.10) rewrite 
formula (1.5) for the transfer matrix function 

wa{x, X)^Im- ijTl{xyS{x)-\A - A7„)-^n(x). (2.11) 

Thus constructed matrix function wa is a gauge transformation of the Dirac 
type system. To show this we shall differentiate Wa using (2.6), (2.7), and 
(2.9). First calculate the derivative of n*5'~^: 

(^n*^-^) = ijii*A*s-^ + iVjivs-^ - ii*s-^im*s-\ (2.12) 

Taking into account that (2.9) yields A*S-^ = S'^A - iS-'^mTS-'^, we 
rewrite (2.12) as 

(^n*5-^) = i]ii*s-^A + [iVj + jirs-^iij - n*5"^n)n*5"\ (2.13) 

Using (2.6), (2.11), and (2.13) one gets 

^WA = -i3[{iVj+jIi*S-^Iij-WS-^Ii)WS-\A-\I^)-^n 

M3n*S-\A - A7„ + XQiA - A7„)-in) (2.14) 
-ijWS-\A - [-i{A- XIn + XQUj - tUjV^ . 

From (2.11) and (2.14) it follows that 

^WA = {tjV + U*S-'U-jU*S-'Uj){wA-In)+Il*S-'U 

+lXjiwA - In) - jIl*S-'llj - {wA - In){lX3 + IjV) 

= (iAj + ijV + WS-^Ii - jIi*S-^Iij)wA - WA{iXj + ijV). 



Finally, rewrite relation above as 



^wa{x, X) ^ i{Xj + jV{x))wA{x, X) - iwA{x, X){Xj + jV{x)), (2.15) 



where 



V 
V* 



V + i{U*S-^Uj - jU*S-^U), 



v^v-2i[Ip Q]U*S-^U 
Thus the following proposition is proved. 



(2.16) 
(2.17) 



Proposition 2.1 Let Dirac system (2.1) and parameter matrices A, 11(0), 
and S'(O) = S'(O)* he given, and let equality (2.5) hold. Then the matrix 
Junction Wa defined by (2.11), where Il{x) and S{x) are obtained via for- 
mulas (2.6) and (2.1), is a gauge transformation (Darboux matrix) of the 
Dirac system and satisfies equation (2.15). The fundamental solution of the 
transformed system -j-u — i{Xj + jV)u is given by the formula 



u{x, A) = wa{x, X)u{x, X)wa{0, A) ^, 
where u is the fundamental solution of the initial system (2.1). 



(2.18) 



Recall that we normalize fundamental solutions u by the condition (1.4). 

The skew-self-adjoint Dirac-typc system can be written as the first rela- 
tion in (2.4), where qi = —ij, qo{x) = —jV{x), and we put 



Ai^A*^ A, Ui = n, = iU*, 



(2.19) 



where only the last equality differs from the last equality in (2.10). After 
substitution of our new qo we define 11 via (2.6): 

U^(x) = AU(x)qi + U(x)qo(x) = -AU(x)j - U(x)jV(x). (2.20) 

We define S by the equality 

= U.q^U* = nj-n*. (2.21) 
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Under these conditions and under the matrix identity at a; = condition 

AS{o) - s{o)A* = m(o)n(o)*, (2.22) 

the identity 

AS(x) - S(x)A* = m(x)U(xy (2.23) 

follows. Similar to the previous Proposition 2.1 our next proposition can be 
proved. 

Proposition 2.2 Let skew-self-adjoint Dirac-type system (2.3) he given, and 
let (2.22) hold. Then the matrix function 

wa{x, X)^Im- ill{,xy S{x)-\A - A4)-^n(x), (2.24) 

where H and S are defined by (2.20) and (2.21), is a gauge transformation 
of (2.3) and satisfies equation 

^wa{x,X) = {iXj+jV{x))wA{x,X)-WA{x,X){iXj+jV{x)). (2.25) 



Here we have 



V 



V 
V* 



V + U*S-^U - jU*S-^Uj. (2.26) 



Proposition 2.2 admits generalization for the case of the system auxiliary to 
nonlinear optics equation 

-^u{x, A) = {tXD - [D, a^)])u{x, A), C = B^B, (2.27) 
where m is an m x m matrix function and 
D^diag{d,,d2,...,drr^}^D*■, 5 = diag {61, 62, ■ ■ ■ , h = ±l- (2.28) 

To present (2.27) as the first equality in (2.4) we should put qi = —iD and 
qoix) — [D,^{x)]. Substitute these expressions for qi and qo into (2.6) to 
define 11, and define S by the equahty Sx = UDBU* for the derivative of S. 
The matrix identity AS{x) — S{x)A* = iIl{x)BU{x)* easily follows after we 
assume AS{Q) - S{Q)A* = in(0)5n(0)*. 



9 



Proposition 2.3 Let system (2.27) be given. Then the matrix function 

wa{x, A) = - iBU{xyS{x)-\A - XQ-^U{x) (2.29) 
is a gauge transformation of (2.27) and satisfies equation 

^M^,^) = {iXD-[Dj{x)])wA(x,X)-WA(x,X){iXD-[D,^(x)]). (2.30) 

Here we have 

^ = BIl*S-^U, C = bIb. (2.31) 

Propositions 2.1-2.3 are particular subcases of a general version of Backlund- 
Darboux transformation for systems depending rationally on the spectral 
parameter A (see Section 3). 

RemEirk 2.4 When A is a scalar (i.e., n — 1), then by (2.6) A is an eigen- 
value and n is an eigenfunction of the dual to (2.4) system Ux = {Xqi -\-qo)u. 
In a general situation we call A a generalized matrix eigenvalue and we call 
n a generalized eigenfunction. 



2.2 N-wave equation and gauge transformation 

Consider N-wave (nonlinear optics) equation 

[D, Ct{x, t)] - [D, U^, t)] = [[D, C{x, t)], [D, C{x, t)]\ (Ct = (2.32) 

on the product Xi x I2 of intervals Xi and X2, where 

e^B^B, (0,0)eXixX2, D ^ dis.g{di, . . . ,dm} ^ D\ (2.33) 

and D and B are given by (2.28). Nonlinear integrable equation (2.32) is the 
compatibility condition of two auxiliary linear systems [21,82] (see also [2] 
for the case N > ?,): 



Ux{x, t, A) = {iXD - [D, ^{x, t)])u{x, t, A), 
ut{x, t, A) = {iXD - [D, ^{x, t)])u{x, t, A). 



(2.34) 
(2.35) 
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Indeed, in view of (2.34) and (2.35) we have 

U:ct = -[D.it]u + {tXD - [D,i]){t\D - [D,i])u, 

and 

ut. = -[5, Qu + {i\D - [D, i\) {iXD - [D, ^])u. 

Thus one can easily see that the compatibihty condition u^t — Utx is equiva- 
lent to (2.32). 

To construct gauge transformation we fix n > and three parameter 
matrices, that is, two nxn matrices A and 5'(0, 0) = 5'(0, 0)*, and an n x m 
matrix n(0, 0) such that 

AS{Q, 0) - 5(0, 0)A* = in(0, 0)5n(0, O)*. (2.36) 

Now, introduce matrix functions n(x, t) and S{x, t) by the equations 

= -iAUD + n[D, {], iit = -lAiiD + n[5, e], (2.37) 

S:, = nDBXr, St = UDBU*. (2.38) 

Quite similar to u^t — Utx one can show that according to (2.32) equations 
(2.37) are compatible, i.e., Hxt — ^tx- By (2.37) and (2.38) it is immediate 
that Sxt — Stx- Proposition 2.3 implies 

Proposition 2.5 Let m x m continuously differentiable matrix function ^ 
(^* — B^B) satisfy N-wave equation (2.32) and let matrix functions 11 and 
S = S* satisfy (2.36)-(2.38). Then in the points of invertibility of S the 
matrix function 

^(x, t) := C{x, t) - BU{x, t)*S(x, t)-^U(x, t). (2.39) 

satisfi.es equation (2.32) and an additional condition ^* = B^B. 

Proof. Equality S,* = B^B is immediate. Now, let u be the fundamental 
solution of (2.34) and (2.35), i.e., let u satisfy equations (2.34) and (2.35), and 
equality -^(0,0, A) = /„• (As ^ satisfies (2.32), the compatibility condition 
for systems (2.34) and (2.35) is fulfilled, and such a matrix function u exists. 
See, for instance, formula (1.6) on p. 168 in [76].) Put 

wa{x, t, X)^Im- iBU{x, tyS{x, t)-\A - XQ-^U{x, t), (2.40) 
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and calculate derivatives of wa with respect to x and t using Proposition 2.3 
in both cases. Then, for u{x, t, A) = wa{x, t, X)u{x, t, A) we have 

Ux{x,t, X) = G{x,t, X)u{x,t, X), Ut{x,t, X) = F{x,t, X)u{x,t, X), (2.41) 

where 

G^iXD-[D,^i F^iXD-[D,l]. (2.42) 

From the continuous differentiability of ^ follows the continuous differentia- 
bility of ^ in the points of invertibility of S. Hence, in view of (2.41) we 
get u^t = or equivalently Gt — F^. + [G, F] = 0. As we have already dis- 
cussed here, the last equality is in its turn equivalent to the N-waye equation 

In the case of the trivial initial solution ^ = we obtain explicit solutions 
of the A^-wave equation. Namely, putting n(0, 0) — [fi f2 ■ ■ ■ fm] and using 
(2.37) we recover 11: 

U{x,t) = exp [-i{diX + dit)A)fi exp [-i{d2X + d2t)A) f2 ... . (2.43) 

Next we recover S, and explicit formulas for solutions of the N-wave equation 
are immediate from (2.39). 



2.3 Nonlinear Schrodinger equation: 
n-modulation solutions 

Auxiliary systems (2.34) and (2.35) are a particular case of the hnear differ- 
ential first order systems 

Ux{x, t, A) = G{x, t, X)u{x, t, A), ut{x, t, A) = F{x, t, X)u{x, t, A). (2.44) 

The compatibihty condition for systems (2.44) is given by the equation 

Gt-F, + [G,F]=0. (2.45) 

When the first and second system in (2.44) takes the form (2.34) and (2.35), 
correspondingly, the compatibility condition (2.45) is equivalent to N-wave 
equation (2.32). 
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In this subsection we consider another example, namely, the well-known 
integrable [85] focusing nonlinear Schrodinger equation (fNLS) 

2vt + i{v^^ + 2vv*v) = 0. (2.46) 

Here v{x,t) is a p x p matrix function. Equation (2.46) is the compatibility 
condition of the auxiliary systems (2.44), where 

G = i\j + jV, F = - iXjV - {V, + jV^)/2) , (2.47) 

and j and V{x, t) are defined in (2.2). Consider again domain Xi x X2, where 
(0, 0) G Xi X I2. Our next proposition is a particular case of the results of 
Section 3 that can be proved similar to Propositions 2.1, 2.5. 

Proposition 2.6 Letpxp matrix function v{x,t) satisfy equation (2.46) and 
be continuously differentiable together with v^. Let n x m matrix function U 
and n X n matrix function S — S* satisfy equations 

AS{0, 0) - 5(0, 0)A* = iU{0, 0)n(0, 0)*, (2.48) 

n, = -iAUj - UjV, = -iA'Uj - AUjV + '-U{V, + jV^), (2.49) 

= UjU*, St = AUjU* + njEM* - tUjVU*, (2.50) 
where A is an n x n parameter matrix. Then the identity 

AS(x, t) - S(x, t)A* = iU{x, t)U(x, t)* (2.51) 

holds in the domain Ti x I2, and the matrix function 



V + 2[ip o]n*5-^n 







(2.52) 



satisfies (2.46) in the points of invertibility of S. Moreover, the transfer ma- 
trix function wa given by (2.24) satisfies equations {wa)x = Gwa — waG, 
iwA)t — Fwa — waF , where G and F are obtained by substitution of V — 
V 1 

^ instead of V into the right-hand sides of the first and second, re- 
spectively, relations in (2.47). 
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In the previous subsection we constructed explicit solutions of the A^-wave 
equation, using trivial initial solution ^ = 0. In the same way the A^^-soliton 
solutions of the fNLS can also be constructed [28,54], putting v = in (2.49), 
(2.50), and (2.52). Now, we shall study a slightly more complicated situation, 
that is, the case of the nontrivial background . 

Example 2.7 Let p — 1, and v — e~'*. One easily checks that 

u{x,t,X) = (exp(-itj72))Co(A)exp ((a; + At)Ci(A)), (2.53) 

where 



Co{X) 



1 1 



(2.54) 

satisfies (2.44), (2-4'^) for the case v — e~'*. Choose for simplicity 

^ = diag{ai,a2, . . . ,a„}, ak^ai {k,l<n), (2.55) 

and put 

Mx,t)^{u{x,t,a^)fky, /ikeC^ /fe^O. (2.56) 
It follows from (2.44), (2-47), and (2.56) that 

irk). = + jvrk, mt = m'^oru + aj^jv^^i - \jv^i^i- (2.57) 

Taking into account that — (i.e., — 0), by (2.55) and (2.57) we see 
that the matrix function 

' ^1 " 

n= ... (2.58) 

satisfies (2.49). Thus H is obtained from (2.56) and (2.58), and then rela- 
tions 

S = {sfej}fcj=i, Skj = ii^ki^*/{ak - a]) (2.59) 

follow from the matrix identity (2.51). In this way, using (2.52) matrix 
function v is constructed explicitly. 
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Fix now integers {rik}, {r2k} such that ~ ^ifc — ^fc ^ k < n), 
where 1^ are integer, and put = iT\k/f2k- Then, by (2.53) and (2.56) the 
dependence ofipl on t can be expressed in terms of functions exp(±it/2) and 
exp{±ilkrikt /r2i^) . Therefore, taking into account (2.52), (2.58), and (2.59) 
we see that v is a periodical in t solution. 

Under somewhat more restrictive than (2.55) conditions 

A = diag{ai,a2,...,a„}, (7(A) e C+, ak^ai{k^l), (2.60) 

matrix function S is invertible, and moreover we have S > 0. To show this 
we rewrite identity (2.51) in the form 

S{A* - XQ'^ -{A- XIn)-^S = i{A - XIn)-^UU*{A* - XQ-'. (2.61) 

As o'{A) G C+, by the theorem on residues one represents S as integrals of 
the right-hand side of (2.61) on contours in C+, and in the limit we get 

1 /"°° 

S^— {A- XIn)-^mr{A* - XQ-'dX, (2.62) 

that is, S >0. To derive the strict inequality suppose Sg — 0, g — {gj}^=i 7^ 
0. Hence we have g*{AS—SA*)g — 0, and so by (2.51) the equality g*TlU*g — 
holds, i.e., IT* (7 = 0. Using again (2.51) one gets SA*g = 0. Now, by 
induction equalities 15(^4*) g = 0, 11* (A*) g = (k > Q) easily follow. As 
by our assumption g ^ 0, there is its entry gr ^ 0. In view of the third 
relation in (2.60) we obtain e(r) := {5r,j}j=i ^ span |J^~q(A*) g, where Srj 
is the Kronecker-symbol. Therefore from Il*[A*) g = (k > 0) it follows 
n*e(r) = 0, that is, 'ipr = 0. The last equality contradicts (2.56), and so 
inequality S > is proved. Taking into account S > we see that the fNLS 
solutions given by (2.52) are well-defined. 
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3 GBDT for system depending rationally on 
spectral parameter and explicit solutions of 
nonlinear equations 

3.1 GBDT for system depending rationally on A 

In this section we consider GBDT for a general case of first order system 
depending rationally on the spectral parameter A: 

r I Ts 

= Gu, G{x, A) = - ( ^ \'qk{x) + X^(A - CsY'^qskix)) , (3.1) 

fc=0 s=l k=l 

X G X, where X is an interval such that G X, and the coefficients qk{x) and 
Qsk{x) arc m x m locally intcgrable matrix functions. 

As before we fix an integer n > 0. Next, we fix five matrices, namely, 
n X n matrices Ak {k = 1, 2) and S{0), and n x m matrices 11^(0) {k = 1, 2). 
It is required that these matrices form an 5'-node, that is, the identity 

Ai^(O)- 5(0)^2 = Hi (0)n2(0)* (3.2) 

holds. Matrix functions 11^(2;) are introduced via the coefficients from G: 

(ni), = ^i^i^fc + Y. $^(^1 - cJn)-'nig3ik, (3.3) 

fc=0 s=l k=l 

r I Ts 

(n;). = - ( E ^^n;^2 + E E <iskmA2 - cjn)-") . (3.4) 

k=0 s=l k=l 

Compare (3.1) with (3.4) to see that Ilg can be viewed as a generalized 
eigenfunction of the system Ux = Gu. 

Matrix function S{x) is introduced via ^'S' by the equality 

r k I Ts k 

= J^J^A^-^nig.n^Ar'- J]X^$^(Ai-cA)^-^-' (3.5) 

k=l j=l s=l k=l j=l 

xUiqskIl*2{A2 - Csln)~^- 
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Equality (3.5) is chosen so that the identity (AiS-SA2^ = (niH*) holds. 
Hence, taking into account (3.2) we have 

AiS{x) - S{x)A2 = Ui{x)U2ixy, xel. (3.6) 

By Theorem 3.1 below, the Darboux matrix for system (3.1) has the form 
(1.5) : 

WAix, A) = /„, - Mxysixy'iAi - xQ-'Mx). (3.7) 

In other words, wa satisfies the equation 

-^wa{x, A) = G{x, X)wa{x, A) — wa{x, X)G{x, A), (3.8) 
where G has the same structure as G: 

G{x, A) = - ( J] X%{x) + E - ^s)-%k{x)) . (3.9) 

A;=0 s=l k=l 

The transformed coefficients qk and Qsk are given by the formulas 

r j 

qk = qk- ^ (^qjYj_k-i - Xj_k-iqj + ^ Xj_iqjYi_k-2^, (3.10) 

j=k+l i=k+2 

Qsk = Qsk (3.11) 

rs j 

~l~ ^ ^ (^sj^s,k—j—l -^s,k—j—lQsj ^ ^ s,i—j —iQsj^s ,k—i—l^ ) 
j=k i=k 

where Xk{x), Yk{x), Xsk{x), and Ysk{x) are expressed in terms of the matrices 
Ak and matrix functions S{x) and ^.^{x): 

Xk = I{lS-^A\lli, Yk = W^A^^S-^Ili, (3.12) 
X,k = TllS-\A^ - Cslnfn^, Y,k = n;(^ - cJnfS-^n^. (3.13) 

Theorem 3.1 [56] Let first order system (3.1) and five matrices S{Q), A^, 
and rife {k — 1,2) he given. Assume that the identity (3.2) holds. Then the 

transfer matrix function w A given by (3.7), where S andHk are determined by 
(3.3)-(3.5), satisfies equation (3.8), where G is determined by the formulas 
(3.9)-(3.13). 
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The proof of Theorem 3.1 for the case of one pole, that is, for G{x,X) — 
Yyk=-r ^''ik is contained in [60] . The case of several poles Cg can be treated 
precisely in the same way. The following formula is essential for the proof 
and is also of independent interest: 

r I rs 

^ A;=0 s=l k=l 

Formula (3.14) means that multiplying from the right we transform 

a generalized eigenfunction of system (3.1) into a generalized eigenfunction 
of the transformed system Ux = Gu. (Compare formula (3.14) with formula 
(3.4).) 

Remark 3.2 It is immediate from (3.10) and (3.12), respectively, thatqr — 
Qr and Xq — Yq. 

Remark 3.3 // cr(Ai) fl (7(^2) = the matrix function S{x) is uniquely 
defined by the matrix identity (3.6). 

Remark 3.4 In the points of the inveriihility of S {and for A ^ o'(^i) U 
a{A2)) the matrix function wa{x. A) is invertible. Indeed, from the realization 
(3.7) and formula (A. 4) it follows that 

wa{x, A)-^ ^Im + U^ixySixy^A"" - XQ-'U,(x), 

where A^ = A\ — HillgS'"^. In view of (3.6) it is immediate that A'^ — 
SA2S~^. Hence, we get 

Wa{x, A)-^ ^Im + - \In)'^S{x)-^Il^{x). (3.15) 

3.2 Explicit solutions of nonlinear equations 

One can apply Theorem 3.1 to construct solutions of nonlinear integrable 
equations and corresponding wave functions similar to the way, in which it 
was done in subsections 2.2 and 2.3. For this purpose we use auxiliary linear 
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systems for integrable nonlinear equation: 

Ux — Gu, Ut = Fu; (3.16) 

r I Vs 

G{x, t,X)^-J2 >^'Qk{x, t)-J2 - ^sr'qskix, t), (3.17) 

k=0 s=l k=l 

R L Rs 

F{x, t,X) = -J2 >^'Qk{x, t)-Y^ - Cs)-'Qsk{x, t), (3.18) 

k=0 s=l k=l 

and zero curvature (compatibility condition) representation (2.45) of the in- 
tegrable nonlinear equation itself. We consider nonlinear equations in the 
domain (x, t) eXi x X2 and assume (0, 0) e Xi x X2. 

Remark 3.5 If G and F are continuously differentiable and (2.45) holds, 
then according to formula (1.6) on p. 168 in [76] there is the m x m solution 
u of (3.16) normalized by the condition u{0, 0, A) — Im- 

Remark 3.6 Usually we shall asume that G and F are continuously differ- 
entiable and that u{x, t, A) is the solution of (3.16) normalized as in Remark 
3.5. 

When we deal with two auxiliary linear systems, the nx n matrix functions 
S and the n x m matrix functions Iljt depend on two variables x and t, 
and the matrix identity (3.2) for parameter matrices A^, llk{0), and S{0) is 
substituted by the identity 

AiS{o, 0) - s{o, o)A2 = ni(o, o)n2(o, o)* (3.19) 

for parameter matrices A^, 11^(0,0), and S'(0,0). Equations (3.3)-(3.5) 
should be completed by the similar equations with respect to derivatives 
in t. Then Theorem 3.1 provides expessions for derivatives {wA{x,t, X)) ^ 
and (wA{x,t, X))^, where 

wa{x, t, x)^im- n2(x, tys{x, t)-\A^ - A4)-^ni(x, t). (3.20) 

Hence, equations 

Ux — Gu, Ut — Fu, u{x,t, X) :— WA{x,t, X)u{x,t, X) (3-21) 
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hold. Finally in a way which is similar to the proof of (3.6), one can show 
that 

AiS{x,t) - S{x,t)A2^Ui{x,t)U2{x,t)*, {x,t)elixl2. (3.22) 

It follows from (3.21) that 

u,t = {Gt + GF)u, ut, = {% + FG)u. (3.23) 

If G and F arc continuously differentiable, then G and F are continuously 
differentiable too. Hence, Utx = u^t and formula (3.23) implies 

{Gt-F, + [G,F])u^O. (3.24) 

By Remark 3.4 and identity (3.22) the matrix function WA{x,t,X) is invert- 
ible, and by Remark 3.5 u is invertible. Thus, u is invertible. Therefore, it 
is immediate from (3.24) that 

Gt-Fx + [G, F] = 0. (3.25) 

For the particular case of the A^-wave equation formula (3.25) was obtained 
in the proof of Proposition 2.5 and was used there to show that ^ satisfies 
the N-wecve equation. Now, we proved the following general theorem. 

Theorem 3.7 Let G and F be continuously differentiable and satisfy (2.4.5). 
Let the identity (3.2) hold and let the matrix functions 11^ and S he given 
by the equations (3.3)-(3.5) and by the analogs of (3.3)-(3.5) with respect to 
the variable t instead of x, namely, by the equations 

R L Rs 

= E ^tniQfe + - c^.4)-'niQ,fe, (3.26) 

fe=0 s=l k=l 

R L Rs 

(n;)t = - ( E ^'^"2^12 + E E Qsk^*2(A2 - CsQ-') . (3.27) 

k=0 s=l k=l 

R k L Rs k 

fc=l j=l s=l k=l j=l 

Then in the points of the invertibility of S the zero curvature equation (3.25), 
where G and F are given by (3.9)-(3.13), holds. 
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Example 3.8 The main chiral field equation for m x m invertible matrix 
function z has the form 

2z^t{x, t) = z^{x, t)z{x, ty^zt{x, t) + Zt{x, t)z{x, ty^z^{x, t), (3.29) 

and is equivalent [49,83] to the compatibility condition (2.45) of the auxiliary 
systems (3.16), where 

G(x,t,X)^-(X-l)-^qii{x,t), qn^z^z-^; (3.30) 
F{x,t,X)^-{X + l)-'Qi,{x,t), Qu^-ZtZ-\ (3.31) 

Let z satisfy (3.29). In view of (3.30) and (3.31) equations (3.3)-(3.5) 
take the form 

(Hi)^ = (A - I^)-^Il^z,z-\ (n;)^ = -z^z-^ni^A^ - Q-\ (3.32) 

Sx = -(A, - Q-^U,z,z-'Ul{A2 - In)-\ (3.33) 
and equations (3.26)-(3.28) take the form 

(Hi)^ = -(^1 + Q-^U,ztZ-\ (n;)^ = Ztz-^UI{A2 + In)-\ (3.34) 

St = (Al + In)'^\l^ZtZ-^\ll{A2 + In)-^. (3.35) 

Now, let matrices Ak, 5'(0, 0), anc? 11^(0, 0) he fixed, assume that (3.19) holds, 
and let matrix functions S andUk satisfy (3.32)-(3.35). By Theorem 3.1 we 

get 

(wa)^ = Gwa - waG, {wA)^ = Fwa - waF. (3.36) 

Assume additionally that det A^. 7^ (/c = 1,2). Then, by (3.15) and (3.20) 
the matrix functions wa^x., t, 0) and wa{x, t, 0)~^ are well-defined in the points 
of iinvertibility of S. Taking into account (3.9) and first equalities in (3.30) 
and (3.31) we have 

G{x,t,X) ^-{X-l)-%i{x,t), F{x,t,X)^-{X + l)-^Qu{x,t). (3.37) 

It follows from (3.36) and (3.37) that 
d 

—wa{x, t, 0) = qii{x, t)wA{x, t, 0) - wa{x, t, 0)qu{x, t), (3.38) 
d 

—WA{x,t,0) = -Quix,t)wAix,t,0) + WA{x,t,0)Qii{x,t). (3.39) 
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Rewrite second relations in (3.30) and (3.31): 

= QuZ, Zt = —QiiZ. 

Put 



(3.40) 



z{x,t) := WA{x,t,Q)z{x,t). (3-41) 

From formulas (3.38)-(3.41) we derive = quz and Zf = —QnZ. As 
WA{x,t,0) and z{x,t) are invertible, so z is invertible, and we get 

'\ Qu = -ZtZ-\ (3.42) 



ZxZ 



Recall that if formulas (3.30) and (3.31) hold, then (2.45) is equivalent to 
(3.29). The only difference between equalities in (3.30), (3.31) and equalities 

in (3.37), (3.42) is "tilde" in the notations. Hence, in view of (3.37) and 
(3.42) formula (3.25) implies that z satisfies main chiral field equation, that 
is, 

2zxt{x, t) = Zx{x, t)z{x, ty^zt{x, t) + zt{x, t)z{x, t)~'^z^{x, t). 

Corollary 3.9 Assume that parameter matrices satisfy identity (3.19) and 
that Act Ak 7^ (A; = 1,2). Let an invertible matrix function z satisfy main 
chiral field equation (3.29) and be two times continuously differentiable. Then 
the matrix function 'z given by (3.41) in the points of invertibility of S also 
satifies main chiral field equation. 

Our next examples deal with the construction of new (local) solutions of 
integrable elliptic Sine- Gordon and sinh- Gordon equations from the initial 
solutions. See, for instance, [10,35] and references therein for some related 
literature and auxiliary systems . 



Example 3.10 Elliptic Sine-Gordon equation 

Vtt + Vxx = sin V {v ^v) 



(3.43) 



is equivalent to the compatibility condition (2.45) of the auxiliary systems 
(3.16), where 



F = -^{XC + Vxj + jJCJ), 







(3.44) 
(3.45) 
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and matrices j and J are defined in (1-7) after putting p — 1. We put also 

= A, A2 = -{A*)-\ Hi = n, n2(0, 0) = A-^U{0, 0) J. (3.46) 

Thus, we have three parameter matrices, that is, nxn matrices A and S{0, 0) 
and annxm matrix 11(0, 0) . We assume that v satisfies ( 3.43), that det A ^ 
and S'(0, 0) = 5'(0, 0)*, and that there is a matrix U such that equalities 

A = uA-^u-\ n(o, 0) = c/n(o, o), s{o, o) = uas{o, o)a*u* (3.47) 

hold. Here A is the matrix with the entries, which are complex conjugate to 
the corresponding entries of A. By (3.46) the identity (3.19), which should 
he satisfied by the parameter matrices, takes the form 

AS{Q, 0)A* + 5(0, 0) = n(0, 0) jn(0, O)*. (3.48) 

Compare (3.17) and (3.18) with (3.44) <^'>T'd (3-45), respectively, to see that 

r = 1, qi = (-V4)C, Qo = {-Vt/4:)j; / = ri = 1, Ci = 0, qn = {i/4:)JCJ; 

R^l, Q,^ (1/4)C, Qo = {vjA)j; L = i?i = 1, Ci = 0, Qn = (1/4) JC J- 
Thus, in view of (3.46) equations (3.3) and (3.26) take the form 

= 1 ( - anc - vtXij + lA-'YiJCJ) , = ^ (Anc + vjij + a~'yijcj) . 

(3.49) 

As C = JC,J o-nd A = UA~^U~^, one can see that both U and z = UU satisfy 
(3.49). According to (3.47) n(0, 0) = I7n(0,0)), and so we derive 

\l{x,t)=Tm{x,t). (3.50) 

Equations (3.4) and (3.27), which define 112, ^^^^ the form 

(n;), = \{iCniA2 + vtjni - iJcmiA-'), (3.51) 
(n;)t ^\{- cn;A2 - v^jw^ - JCJn^^-^). (3.52) 

^5 C = C* and A2 = -{A*)-\ it follows from (3.49) that 

z{x,t) = A-^U{x,t)J 
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satisfies equations (3.51) and (3.52) for 112. Moreover, we have 112(0,0) = 
74~^n(0, 0) J = 2;(0, 0). In other words we have 

Ii2{x,t) = A-^Ii{x,t)J, (3.53) 

and identity (3.22) takes the form 

AS{x, t)A* + s{x, t) = n(x, t)m{x, ty. 

By (3.5), (3.28), (3.53), and by the second equality in (3.46) the relations 

= ^-[A-^UJCU* - UCJU*{A*)-^), St = ^{A-^UJCU* + UCJU*{A*)-^) 

(3.54) 

hold. Formulas (3.47), (3.50), and (3.54) imply = UAS^A*U*, St = 
UAStA*U* and finally 

S=UASA*U*. (3.55) 

It follows from (3.12), (3.50), (3.53), and (3.55) that 



x_i = m*{A*)-^s-^A-^Ti = m*u*{uASA*u*)-^un = j(n*,s-in). 

(3.56) 

According to (3.12), (3.15), and (3.20) we have 

Z{x, t) := wa{x, t, 0) = 72 - ^-1, Z{x, t)-^ = wa{x, t, 0)"^ ^ I2 + Y_i. 

(3.57) 

Moreover, in view of (3.12), (3.53), and equality U*S-^U = (n*,S-^n)* we 
get 



y_i = -ju*s-^u, -u*s-'u 



Using (3.56)-(3.58) we derive 



a b 
b d 



, a = a, d = d. (3.58) 



1 + 6 d 
a 1 + b 



1 + b d 
a 1 + 6 



If 1 + b ^ 0, formula (3.59) implies 

a^d = Q, |1 + 6| = 1, Z = diag{l + 6, 1 + 6}. 



(3.59) 



(3.60) 
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Put 



f + 2arg(l + 6), u = Z ^u, G = UxU 



UtU 



(3.61) 



where u{x, t, A) = wa{x, t, X)u{x, t, A). In a way, which is similar to the proof 
of (3.25) in Theorem 3.7, we derive from (3.61) that 



Gt-F, + [G,F]=^. 



(3.62) 



Moreover, one can see that G and F have the form (3.44) '^'f^d (3.45), respec- 
tively, after one substitutes v instead of v into the right-hand sides of (3.44) 
and (3.45). Therefore, formula (3.62) implies thatv satisfies (3.43). 

Corollary 3.11 Let an integer n > and matrices A (detA ^ 0), 11(0,0), 
and 5'(0, 0) = 5(0,0)* be fixed and satisfy conditions (3.47) and (3.48). Let 
V satisfy elliptic Sine- Gordon equation and be two times continuously diff'er- 
entiable. Then in the points, where det S and 1 + 6 7^ 0, the function v 
given by (3.61) satisfies elliptic Sine-Gordon equation too. 

If a (A) n a((— A*)^^) = 0, then the last equality in (3.47) follows from 
(3.48) and from the first two equalities in (3.47). 

Remark 3.12 Using considerations from (3.57) one easily shows that a gen- 
eral equality [Im — -^-i)(/m + ^-1) = Im is true. 



Elliptic sinh- Gordon equation 

Vtt + = sinh V 



[V = V) 



(3.63) 



is equivalent to the compatibility condition (2.45) of the auxiliary systems 
(3.16), where 



F = -i(-^AC + ^^.J + ^C)■ 





,v/2 



-v/2 



(3.64) 
(3.65) 



Put 



= A A2^-{A*)-\ ni = n, n2(o,o) = A-in(o,o). (3.66) 
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Here we assume that 

det A^O, S{Q, 0) = S{Q, 0)*, (3.67) 
and that there is a matrix U such that 

A = UA-^U-\ n(0, 0) = t/n(0, O) J, :S(0, O) = UAS{Q, Q)A*U*. (3.68) 
Now, the identity (3.19) takes the form 

A5(0, 0)A* + 5(0, 0) = n(0, 0)n(0, O)*. (3.69) 
Taking into account (3.64) and (3.65) introduce 11 by the equations 

n, = ^ (^nC - ivtHj + A-^nc) , ^t^\{- iAHC + Ivjlj + iA-^nc) . 

(3.70) 

It is easy to see that 112 = ^"^H, and so formulas (3.5) and (3.28) take the 
form 

= -(ncn*(A*)-i + A-^nc*n*), St = -{a-^ucu* -ucu^a*)-^). 

^ ^ (3.71) 

The matrix Z — I2 — X^i is again a diagonal matrix, and we have 

Z^h- Tl*{A*)-^S-^A-^Ii = diag{Zn, Z^Y}, Zn = Z^^. (3.72) 

(compare with formula (3.60), the proof of (3.72) is similar). The following 
corollary is proved in a quite similar way to the Corollary 3.11. 

Corollary 3.13 Let an integer n > and matrices A, 11(0, 0), and 5(0, 0) he 
fixed and satisfy conditions (3.67)-(3.69). Let v satisfy elliptic sinh-Gordon 
equation and he two times continuously differentiahle. Then in the points, 
where det 5 7^ and Zu ^ 0, the function 

?? = t; + 21n|Zii| (3.73) 

satisfies elliptic Sine-Gordon equation too. Here Zu is given hy (3.70)- 
(3.72). 
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4 GBDT for radial Dirac equation 



Radial Dirac equation has the form 

d K 



or equivalently 



{ — 'i<^2-, — I — c"i + Vix)\u — Xu (x > 0), 
\ ax X / 

(^-^ + Xqi + qo{x)^u{x,X) ^0 {x > 0), 



(4.1) 



(4.2) 



where k is integer, cTj are Pauh matrices, 



"01" 




" -i " 




" 1 " 


10 




_ i 


, 0-3 = 


-1 _ 



V{x) = Ve{x)l2 + Va{x)ai + f s(x)(73. 



(4.3) 
(4.4) 



Ve, Va, and Vg are real-valued functions, which are locally integrable on [0, /) 
(here and further I e R+), 



qi = -J, J := {-ia2) 



1 
-1 



(4.5) 



go(a;) = f*(a;)cr3 + J(ve{x)l2 + t's(.T)a3). t'*(ic) = - + fa(a;). (4.6) 

a; 

Here Ve, Va, and represent the electrostatic potential, the anomalous mag- 
netic moment, and the sum of the mass and the scalar potential, respectively. 

If we put p = 1 in (1.7) we have cri = J and = j. Recall that Dirac- 
type system of the form (2.1), (2.2) was treated in Subsection 2.1. The radial 
Dirac equation differs from the Dirac-type system. Its structure (as well as 
the structure of V in this section) is somewhat different and it usually has 
singularity ai x — 0, which is of interest from the physical point of view. We 
consider equation (4.1) independently from the results in the Subsection 2.1. 

The double commutation method was applied to (4.1) in an interesting 
paper by G. Teschl [79]. By this method S. Albeverio, R. Hryniv, and Ya. 
Mykytyuk [4] proved that 1 is added to x, when an eigenvalue is removed 
and that 1 is subtracted from x, when an eigenvalue is inserted. (In this 
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section we actively use some of the results from [4].) We apply GBDT to the 
radial Dirac equation (4.1). In particular, we construct explicitly potentials 
and fundamental solutions for the equation (4.1) with x > starting from 
the trivial equation (i.e., from qo = 0). The case formally corresponds to the 
removal of the eigenvalues, the iterated double commutation formulas for the 
insertion of the eigenvalues are given in [79] . 

Fundamental solution u in this section is a non-degenerate 2x2 solution 
of (4.1), we do not require u to be normalized at a; = 0. 



4.1 Main result 



The following procedure to construct explicit solutions of the radial Dirac 
equation is an immediate corollary of Theorem 4.9 from Subsection 4.2. 

Theorem 4.1 To construct a class of explicit solutions of equation (4-1) 
with some fixed integer n, fix an integer m > and m x m matrices Ai and 
»Si > 0. Fix also a >c x x lower triangular matrix A2 and an {m + x) x 2 
matrix n(0) = col[\l/i(0) \E'2(0)], where x = \k\, col means column and ^1 
(^2) is an m X 2 upper (x x 2 lower) block ofU. It is required that 

AiSi - SiAl = ^i(O) J^i(O)*, ^2(0) J^2(0)* = 0. (4.7) 

Moreover, for h := [1 . . . 0]^2 ''^e assume that h{0) = c[0 1], when 
K is positive and odd or negative and even and that h{0) = c[l 0], when k 
is positive and even or negative and odd (c 7^ 0). 
Introduce matrix A 



A := 



R 




A2 



i? = *2(o)j*i(o)*5r\ 



and vectors 9i, 62 ^ £m+>c 

[9, ^2]:=n(0)X, 

Now, put 



K:= — 
V2 



U{x) := [e^^^^i e 



S{x) 



"5i 


" 




_ 








itA 



9i9le 



-itA* + e-''^929;e''^' 



(4.8) 

(4.9) 

(4.10) 
dt. (4.11) 
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Then, if S{x) > for x > 0, the potential 

%{x) := JX{x)J* - X{x), X = {Xi,}].^^ := n^^'^n (4.12) 
admits representation 

%{x) = -(73 + T(x) = -(73 + Va{x)(T2, - Vs{x)(Ti, (4-13) 
X X 

where Va, and Vg are bounded in the neighborhood of zero, and 

Va^X22-Xn--: Vs^X,2 + X2,. (4.14) 

X 

The fundamental solution u{x, A) of system 

+ Ml + qo{^))u{x, \) = {x>Q), (4.15) 
is given by the formula 

u{x,X) = wa{x, X)K expi^—iXxa^), (4.16) 

where 

wa(X) ^h- m*S-\A - A7„)"'n. (4.17) 

By the neighbourhood of zero we mean the neighbourhood of the form (0, e) 
or [0, e). Note that if S{x) > for all x e (0, e), then S{x) > for all 
X e R+, the formulas in Theorem 4.1 are well-defined on 1R+, and is 
infinitely differentiable on M+. 

4.2 Superposition of Darboux transformations 

In this subsection we apply GBDT to the radial Dirac equation. We factorize 
also the Darboux matrix, so that GBDT can be cosidered as a superposition 
of two other GBDTs. Finally, we formulate a general Theorem 4.9. 
In view of (4.2), (4.5), and (4.6) we get 

r = 1, = -J; ql = -JqkJ-\ A; = 0, 1, J"^ = J* = -J. (4.18) 
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Here we fix n > and parameter matrices Ai — A, Il{xo), and S{xo) — 
S{xo)*. Formula (3.3) for IIi = 11 takes the form 

U^ = AUq. + Uqo. (4.19) 

Putting 

A2 = A\ U* = m* (4.20) 

and taking into account equation (4.19) and the third equahties in (4.18), we 
sec that equation (3.4) for is satisfied. By (4.19) and (4.20) formula (3.5) 
takes the form 

= nn*. 

Thus, we have S{x) — S{x)*. It is required that 

AS{xo) - S{xo)A* = n(a;o) Jn(a;o)*, 

< xq < I. The identity 

AS{x) - S{x)A* = U{x)JU{xy 

follows from (4.19), (4.21), and (4.22). 

Remcirk 4.2 IfU can be continuously extended to n(0); then S can also be 
continuously extended to <S'(0). So, in that case we consider 11 and S defined 
on [0, I), and for (4-19) to be true it suffices that (4-22) holds for Xq — 0. 

The next Corollary of Theorem 3.1 and formula (3.14) is immediate 

Corollary 4.3 Let relations (4-19), (4-21), and (4.22) hold. Then, in the 
points of the invertibility of S{x) [x > 0), we have 

(^S{x)-^U{xj^ (x) = A*S{x)-'^U{x)qi{x) + S{x)-'^U{x)qoix), (4.24) 

where 

go go + JXJ* -X, X := J-^Xq = H^^'^n. (4.25) 
Moreover, the matrix function wa given by (4-iV satisfies the equation 

—WAix, A) = G(x, X)wa(x, A) — waIx, X)G(x, A), (4-26) 
ax 
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(4.21) 
(4.22) 

(4.23) 



G{x, A) = -Agi - go(r), G{x, A) = -Agi - qoix). (4.27) 
In other words, we have 

-^u(x, A) = G(x, X)u(x, A), (4.28) 
ax 

where 

u{x, A) = wa{x, X)u{x, A). (4.29) 

Recall that by Remark 3.2 qi — qi- This explains why we have the coefficient 
— gi in the expression for G in (4.27). Equation (4.28) follows from (4.2) and 
(4.26). This equation is a GBDT transformation of the equation (4.2) and 
has the same structure. Namely, to get go instead of go one substitutes into 
(4.6) -u* and Vg instead of and Vs, respectively. That is, by (4.6) and (4.25) 
we have 

go(a;) = Ve{x)J + v^{x)a3 - Vs{x)ai, (4.30) 

where Ve{x) = Ve{x), 

v*{x) ^ v^{x) + X22ix) - Xii{x) , Vs{x) ^Vs{x)+Xi2(x) + X2i{x), (4.31) 
and Xkj are the entries of X. 

Remcirk 4.4 When n — 1 and A & W our transformation coincides with 
the double commutation transformation for Dirac equations treated in [4, 79]. 
When ^4 ^ R our transformation somewhat differs from the transformation 

in [4], because the transformation in [4] uses the equivalent (forn = 1) o/IT 
and of the transposition ofH, whereas we use here 11 and H* , so that Va and 
Vs are real valued. 

We will be interested in the transformation of the equation (4.2) with k, — 
into equation with integer nonzero k. 



Example 4.5 Let n = 1, qo = 0, S{Q) = 0. Notice that we have 

k* = k-\ (4.32) 



qi = -J = ikaak*, K := ^ 

v2 



1 1 

—i i 
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and put (nK)(0) = [1 a\. Taking into account — Q and (4-32), we 
rewrite (4-19) and (4-21) in the form 



{UK):, = lAIiKai, S:, = {UK){UKy 
It is immediate from (4-33) that 



(4.33) 



(4.34) 



Recall that S{0) — 0. Thus, when A ^ A, we require additionally \a\ — 1 so 
that (4-23) holds at x = 0. Hence, formula (4-23) is true for all x > (see 
Remark 4.2). Using (4.23), (4-33) and (4-34), we get 



S{x) = i{A - A)-' (^e'^(^-^) - e^^(^-^)) for A ^ A, 
S{x) = (1 + \af)x for AeR. 

By (4.31), (4.34) and (4-35) we have for A ^~A the equality 



(4.35) 



)"'. (4.36) 



In a similar way we have 



VAX 



)^-(ae-^'^^ + ae^'^^ 



) (^(1 + |q;|^)x) ^ for Ae 



From (4-36) and (4-37) we derive that 



VAX) 



X 



a-\- a 



(4.37) 



(4.38) 



where Va is continuous on [0, 00). When a = ±1, we obtain k. = ^1, i.e., k 
is integer. 

To study GBDT we will need, in particular, to split it into a superposition of 
several transformations, and our next result is dedicated to this procedure. 
Let y4 be a block lower triangular matrix: 



All 
A21 A22 



S 



Sii Si2 

S21 S22 



(4.39) 
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where An and Sn are rii x rii matrices, A22 and 5*22 are ^2 x ^2 matrices, 
n — rii + n2- In our further considerations wc fix some value of x and 
omit temporarily the variable x in the notations. Assume det 5* 7^ and 
det 5*11 7^ 0, and denote by T22 the n2 x n2 right lower block of T = S~^. The 
invertibility of T22 follows from the invertibility of 5*11 and S. One can check 
directly that 



^11 ''r ^11 S12T22 S21S11 —Si^Si2T22 

—T22S21SII T22 



T1 — 1 o C c^i c 

22 ~ '-'22 — >-'21'-'ll '-'12- 

(4.40) 

As A is a block lower triangular matrix and det 5 7^ 0, det^n 7^ 0, so wa 
admits factorisation [75, 76] 



WAiX) = W2{X)Wi{X), 



(4.41) 



Wi{X) = /2 - JU*P*S^,\A,^ - A7„J-iPin, (4.42) 

W2{x) = / - m*s-^p;{A22 - xin,)-^T22^P2S-^u, 

where 

Pi = [ ] , P2 = [ In,]. (4.43) 

One can easily see that by (4.39) we have PiA = A^Pi and AP2 = P2*^22- 
Hence, the identity 

AnSn - SiiAl, = njirl m := PiH (4.44) 

follows from (4.23). Now, compare (4.17) and the first relation in (4.42) to 
see that Wi — wah- That is, we get Wi after substitution of An instead of 
A, of Sii instead of S and of tti instead of 11 into (4.17). Identity (4.44) is 
the equivalent of (4.23) written for the new transfer matrix function w^n- 

Moreover, from (4.23) it follows that TA - A*T = TUJU*T. Therefore, 
we derive ^22^22 - ^22^22 = P2TUJWTP* or, equivalently 

^2T22'-T22^^;2=7r2J7r2*, 7r2 = T2-/P25-^n. (4.45) 

Compare (4.17) and (4.42) to see that W2 = wa,,, where is substituted 
into (4.17) instead of 5" and 7r2 is substituted instead of 11. 
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Remark 4.6 Corollary 4.3 defines by formulae (4.25) and (4.29) GBDT 

transformation of the coefficient Qq and solutions u of the radial Dirac equa- 
tion (4-2). The next proposition shows that the GBDT generated by the 
parameter matrices A, S{xo) and n(xo) can be treated as a superposition 
of two transformations, i.e., GBDTs generated by An, S'ii(a;o), 7ri(a;o) and 
^22, T22{xo)~^, 7ri(a;o), respectively. In particular, our next proposition shows 
that TTi and Sn satisfy analogs of (4-19) and (4.21), respectively. The same 
is true for T22 and 1:2- According to Corollary 4-3 it means that A) 
and W2{x, A) are Darboux matrices. 

Proposition 4.7 Let relations (4-19), (4-21), and (4-22) be valid and let A 
be a block lower triangular matrix. Then, in the points of the invertibility of 
S{x) and 5*11(0;) {x > 0), we have 

{Sll}^ = TTiTT*, (7ri)a; = AuTTiQi + TTiQa. (4.46) 
(^22^)^ = TTaTr^, {■K2)x = ^227r2gi + 772^0, (4-47) 

where tti and 1:2 are given by (4-44) and (4-45), respectively, and 

go = go + JXJ* -X, X = 7r*5f/7ri. (4.48) 

Proof. Multiply from the left both sides of (4.19) by Pi and use PiA — 
AuPi to get the second relation in (4.46). Multiply by Pi from the left and 
by P* from the right both sides of (4.21) to get the first relation in (4.46). 
In view of (4.21) and of the definition of tt2 in (4.45), we have also the first 
relation in (4.47): 

d rp-i _ T-^f ^ T \t-^ — T-^p f ^ t\p*t-^ 

^-'22 — ~^22 y^^22jJ-22 ~ ~ 22 -^2^^^-' J -^2 22 

= T-2^P2S-^mVS-^P;T-; = 7r27r;. (4.49) 

Now, use (4.24), (4.49), and equality P2A* = AI2P2 to differentiate TTa: 

{'K2)x = 7r27r2*P25-'n + Ti^^^ P2{A* S-^Tlqi + S-^Tlq^) 

= 7r27r2*P2>S-^n + T2^^P2>S-^ngo + T2^'A;2^2-5-'ngi. (4.50) 
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Taking into account both relations in (4.45), both relations in (4.25), and 
qi — J*, rewrite (4.50) in the form 



(tTs). = 7l27i;P2S-'U + 712^0 + A22T^f P2S-'Uqi 

-n^-hlP^S-^Hqi = ^22 772^1 + 7i2{qo + Jn*S-'UJ* 



(4.51) 



Finally, note that according to the definition of and formula (4.40) we get 



n* - TT^Pa 
It follows that 



n* 4 + 



'12 



-I 



n2 



<Sii^[Sn Su]. (4.52) 



(n* - n;P2)S-^U = 7ri*5fiVi (4.53) 
Substitute (4.53) into (4.51) to derive 

{7T2)x = A22-K2qi + 7^2 (^0 + JTTi^fiVi J* - TTi^fiVi) . (4.54) 

Formulas (4.48) and (4.54) imply the second relation in (4.47). ■ 

Proposition 4.7 deals with the S'-nodes, which appear in the process of fac- 
torization (4.41) of the transfer matrix function wa- An inverse in a certain 
sense result is given in the next proposition. 

Proposition 4.8 Letwj^^. (j = 1,2) he transfer matrix functions of the form 
(4-i7) corresponding to the S -nodes Aj, Sj and where Ai and Si are 
m X m matrices, A2 and S2 are x x x matrices, \E'i is an m x 2 matrix, ^^2 
is a K X 2 matrix. Then we have w_42(A)w^^(A) = u'^(A), where Wa is the 
transfer matrix function corresponding to the S-node of the form 

R = ^fJ^flSi\ 
(4.55) 

\E'jJ\&* {p = 1,2) and representation (4.55), 
where R = ^2 J*I5r\ it follows that AS - SA* = UJU*. Now, from (4.41), 



' Ai " 




'Si " 




" ^1 " 




R A2 _ 




^2 . 


^2 . 



Proof. From AjSj - SjA* 
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(4.42) and (4.55) we have ^^(A) = W2{X)wA,iX). By (4.40), (4.42), (4.45) 
and (4.55) we obtain — S2, 1^2 — ^2 and, finally, W2 — wj\^^. ■ 

Assume for simplicity that go hi (4.2) is bounded in the neighbourhood of 
zero. Using Corollary 4.3 and Proposition 4.7 we shall prove in the next 
section the following theorem. 

Theorem 4.9 Let the initial system (4-2), where the coefficient qo{x) is 
bounded in the neighbourhood of zero, be given. Let the parameter matri- 
ces A, S{0), and n(0) have the block form 



Ai 
R A2 



S{0) 



5i(0) 




n(o) 



^i(O) 

^2(0) 



(4.56) 



where Ai and iSi(O) are m xm matrices, A and S{0) are (m + x) x (m + x) 
matrices {x > 0), ^i(O) is an m x 2 matrix, \E'2(0) is a x x 2 matrix, A2 
is a lower triangular x x x matrix, and R = ^2(0) J^i(0)*5i(0)~^. Let the 
relations 

AiSiiO) - Si{0)Al = *i(0) J*i(0)*, ^2(0) J*2(0)* = 0, 5i(0) > 

(4.57) 

hold. Introduce Il{x) andS{x) by (4.19) and (4-21), respectively, and assume 
S{x) > for X > 0. Put h := [1 . . . 0]^2 O'^d suppose 

/i(0) = c[l 0] or /i(0)=c[0 1] (c^O). 

Then the transformed system 



d 

dx 



+ Agi + qo{x)]u{x, A) = 0, 



(4.58) 



where go is given by (4-25), has the fundamental solution u — wau. Moreover, 
the coefficient % admits representation 



K 



qo{x) = - + T{x), 



X 



(4.59) 



where T is bounded in the neighbourhood of zero. Here, n = x, if x is odd 
and h{0) = c[0 1] or if x is even and h{0) = c[l 0]. We have k — —x, if 
X is even and h{0) = c[0 1] or if x is odd and h{0) = c[l 0]. 
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4.3 Proof of Theorem 4.9 



When n(a;) is squarely integrable in the neighbourhood of 0, using (4.21) we 
get S{x) = S{0) + J;^Il{t)U{tydt. When S{0) = 0, we have 

S(x) = [ Il(t)U(t)*dt. (4.60) 
Jo 

Lemma 4.10 Letn^l, Ae C, S{0) = 0, and n(0) = ai[l a2\K*, where 
«! 7^ and a2 — ±1. Assume that the potential Qq of the initial system (4-2) 

is 

(i) hounded on (0, e) for some £ > 
or 

(a) has a hounded limit, when x tends to +0. 

Then for the GBDT tranformation of Qq defined hy (4-25) we have 

qo(x) = T±(x) T -as, (4.61) 

where Tj^{x) and T_(,t) are hounded in the neighhourhood of zero, if (i) is 
fulfilled, and have a bounded limit, when x tends to +0, if (ii) is true. 

Proof. Note that under conditions of the lemma we have AS (0) — S {Qi)A* — 
0. According to (4.32) we have also n(0) jn(0)* = 0. So, the identity (4.23) 
holds at x = 0, and therefore it is satisfied at all a; > 0. Hence, the matrix 
function % given by (4.25) is the GBDT transformation of go- 
By the definition of K in (4.32) we get 

n(0) = c[l 0] for a2 = 1, n(0) = c[0 1] for as = -1 (c ^ 0). (4.62) 

By (4.19) one can see that 

n(,r) - n(0) = xf{x), (4.63) 

where / is bounded on [0, e), if go is bounded, and / has a bounded limit for 
X — > +0, if go has a bounded hmit. From (4.60), (4.62), and (4.63) it follows 
that 

S{x) = \c\'^x + x^fi{x), (4.64) 
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where /i is also bounded on [0, e), if qo is bounded, and fi has a bounded 
hmit for x — > +0, if go has a bounded hmit. Taking into account (4.25) and 
(4.62)-(4.64) we obtain the statement of the lemma. ■ 

Recall the notations from Proposition 4.7 and consider the case n > l,ni — 1, 
n2 — n — 1, so that tti is the first block (i.e. the first row) of 11, and 7r2 is 
given by (4.45). 

Lemma 4.11 Let the initial system (4-2) be given and let the potential qo be 
bounded in the neighbourhood of zero. Suppose that A is a lower triangular 
nx n {n> 1) matrix, that relations (4-19) and (4-21) hold and that 

5(0) = 0, n(0) jn(0)* = O, S{x) > O for x > O. (4.65) 

Then, we have 

$o(x) = t+(x) - -(73 for 7ri(0) = c[l 0]; (4.66) 

X 

qo{x) = T_{x) + -a^ for 7ri(0) = c[0 1], (4.67) 

X 

where c e C\{0}, % is given by (4-48) and T± is bounded in the neighbour- 
hood of zero. Moreover, the matrix function 7r2{x) is continuous at zero and 
the relations 

lim (T22^) (x) = 0, (T^2^){x)>Q for x > 0, 7r2(0)J7r2(0)* = 

(4.68) 

hold for both cases (4-66) and (4-67). The matrix functions T^^ and Ti2 
satisfy (4.47). 

Proof. First, compare Proposition 4.7 and Lemma 4.10 to see that % from 
(4.48) coincides with go from (4.61). Thus, (4.66) and (4.67) are immedate 
from Lemma 4.10. 

As S{x) > 0, so we get T{x) > 0, T22{x) > and, finally, {T^^){x) > 0. 
Similar to the case treated in Lemma 4.10 we have AS{0) — S{0)A* = 
and n(0)Jn(0)* = 0, i.e., (4.23) holds at x = 0. Hence, the conditions of 
Proposition 4.7 hold. By Proposition 4.7 722^ '^2 satisfy (4.47). 
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Recall that S{0) = 0. Therefore, using formulae (4.60) and (4.64) we 
derive 

^22(0) = 0, lim S2i{x)Su{x)-^Si2{x) = 0. (4.69) 

The second equahty in (4.40) and formula (4.69) imply the first relation in 
(4.68). Prom the definition of 7r2 in (4.45) and from the first relation in (4.40) 
it follows that 

7r2 = [--S2i5iY 4-i]n. (4.70) 

Now, use again (4.60) and (4.64) to see that 7r2 has a limit, when x tends 
to +0. Therefore, the third relation in (4.68) is immediate from the first 
relations in (4.45) and (4.68). ■ 

Next, consider the case of the initial system with k 7^ 0, where V e L2x2(0, £) , 
that is, the entries of V belong to in the neighbourhood of zero. 

Lemma 4.12 Let the initial Dirac system (4-2) he given, where 

qo{x) = i(72 (V{x) + -cTi) = ia2V{x) + -(73, k 7^ 0, (4.71) 

and V G L2x2(0,e). Suppose that A is a lower triangular n x n {n > 1) 
matrix, that relations (4-19) and (4-21) hold, that S and IT are continuous 
at zero and that (4-65) is true. Then, putting Ui = 1, we get 

qo(x) — (73 + T+(x) for K > 0, qo{x) — (73 + T_(x) for k < 0, 

X X 

(4.72) 

where go is given by (4-48) and T± G i^2x2 neighbourhood of zero. 

Moreover, the matrix function Ti2{x) is continuous at zero and the relations 
(4-68) hold. The matrix functions T22 and 712 satisfy (4-4'^)- 

To prove this lemma we shall need Lemma Al from [4] : 

Lemma 4.13 Let Qq be given by (4-'^l)- Assume k G Z\0, A G C, and 
V G ^^2x2 (1 ^ P < 00) '^n the neighbourhood of zero. Then, there are two 
types of nontrivial solutions y — col[yi ^2] of (4.2) . Namely, as x ^ 
+0, either the limit limx^^o x~'^y{x, X) ^ exists and x~'^~^yi{x, X) G L^ 
in the neighbourhood of zero or the limit \imx^+Qx'^y{x, X) ^ exists and 
x'^~^y2{x, X) & LP in the neighbourhood of zero. 



39 



A similar result is true for the bounded V. 

Lemma 4.14 Assume that n G Z\0, A G C, and V G L'^-^^y that is, V is 
bounded in the neighbourhood of zero. Then, the entries of the bounded in 
the neighbourhood of zero solution y of (4-2) have the property: 

If K > 0, then x'''^''^yi{x, X) G L°° in the neighbourhood of zero, and if 
K <0, then x'^'^y2{x, X) G L°° in the neighbourhood of zero. 

Proof. Lemma 4.14 easily follows from the proof of Lemma Al [4]. Con- 
sider, for instance, the case k, > 0. When V — 0, equation (4.2) has a solu- 
tion (f)Q — col[0? 02]) such that the functions x~'^~^(f)'({x, X) and x'''^(f)2{x, A) 
are bounded in the neighbourhood of zero, and (4.2) has another solution 

ijjQ = col['?/'° 1IJ2], such that there is a limit \imx^^Qx'^ijjo{x, X) 7^ 0. By the 
considerations of [4], the bounded solution y of the Dirac equation (4.2) with 
a bounded V satisfies the integral equation 

{I-C)y = ^, C := M^, A) HcPUt, A) 0°(t, A)]y(t) • dt, (4.73) 

Jo 

7(x, A) := 0o(a;, A) (^1 - y [^?(t,A) ijj^it, X)]V{t)u{t, X)dt^, Q < x < e. 

° (4.74) 
It follows from the properties of 4>o, Lemma 4.13 and definition (4.74) that 

x-'^-Si(2^,A) gL~, x-'^72(x,A) gL°° (7 =[71 72],0<x<£). 

(4.75) 

Using (4.73) and (4.75), by induction we get for some C > and e{C) > 
that 

II (^£^7^ {x)\\ < C^■+^a;^■+^ J > 0, 0<x<s, (4.76) 

where || ■ || is the usual vector norm. Asy — Yl^o^''^^ follows from (4.75) 
and (4.76) that x~'^~^yi{x, A) G L°°. For the case k, < the proof is similar. 
■ 

Notice that from S{x) — U{t)U{t)*dt > one easily gets 7ri(x) ^ in 
any neighbourhood of zero. By (4.19) and by the third equality in (4.18) we 
see that .Ptti satisfies (4.2), where A = A^^. Therefore, the next corollary 
immediately follows from Lemmas 4.13 and 4.14. 
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Corollary 4.15 Let the conditions of Lemma 4-1^ be fulfilled, where 1 < 
p < oo. Then for the entries hi and /i2 o/tti = [hi h2\, for some e > and 
for some c e C\0 we have 



lim x-''hi(x) = c ^ 0, x~''-^h2(x) e L^fO, e) if k> 0, (4.77) 

lim x''h2{x) = c 7^ 0, x^-^hi{x) E LP{0,e), if k<0. (4.78) 

Proof of Lemma 4.12. In a similar to the proof of Lemma 4.11 way one 
shows that the conditions of Proposition 4.7 are fulfilled and so (4.47) is true. 

From S{x) > it follows that {T.^^^{x) > for a; > 0, i.e., the second 
relation in (4.68) is valid. Taking into account the equality S{0) = 0, formula 
(4.46) and Corollary 4.15 we see that 

-5ii(x) = y"''7ri(i)7ri(i)*di= (2|A;| + l)-^|c|Vl'=l+i + o(a;2|'=l+i), x ^ +0. 

° (4.79) 
Moreover, in view of (4.60) and Corollary 4.15 we obtain 

S^i{x) = P2ll{t)7ii{tydt = {\k\ + l)-^a;"^l+iP2n(0)7ri(0)* + o(a;l*^l+i^ 

° (4.80) 
Using the second relation in (4.40) and formulae (4.79) and (4.80), we get 
the first relation in (4.68). From the definition of tti in (4.44) and formula 
(4.70) it follows: 

7r2 = -,52i5nVi + P2II. (4.81) 

By (4.77)-(4.81) the matrix function 7:2 is continuous at zero. Hence, the 
third relation in (4.68) follows from the first relations in (4.45) and (4.68). 

It remains to prove (4.72). For this purpose we shall follow a nice scheme 
from [4]. According to (4.77) and (4.79) we have 

JnlSii'iriJ* - itISii'tti = -^(73 + Ti = -^^3 + T2 for k>0, (4.82) 

where Ti, T2 G I'2x2- (Iii ^^e proof of this lemma functions (1 < j < 8) 
are considered in the neighbourhood of zero.) According to (4.78) and (4.79) 
we have 

JttISii'ttiJ* - TTlSii'm = ^^3 + T3 = ^^3 + T4 for < 0, (4.83) 

'-'11 --"ii 
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where e L|x2 (j =3,4). In view of (4.46), rewrite ninlSi-f in the form 

(4.84) 

From (4.48), (4.71) and (4.82)-(4.84) it follows that 

go(a;) = -^(73 - + T5(a;) for k > 0, (4.85) 

go(2:) = ^—^(73 + ^y^a3 + T6(a;) for «: < 0, (4.86) 
where Tj e L^^g (j = 5, 6). Put 

f(x) := a;-2|*^l7ri(a;)7ri(a;)* = a;-'l'=l7ri(a;) JJ*7ri(a;)*. (4.87) 

Recall that Jvr* satisfies Dirac equation (4.2) with go of the form (4.71), 
where the entries of V belong L^. Hence, differentiating (4.87) and taking 
into account Corollary 4.15 we get 

^/(^) -^M^)M^r + x-'\'\m(x)j(^ - ^a,)rm(xr + T,(x), 

where T7 e i/. Apply again Corollary 4.15 to the right-hand side of (4.88) 
to obtain 

^/(^) = - + ^^^i(^)^3vri(a;)* + Tr{x) = Ts{x), 

(4.89) 

where Tg G L^. Next, use (4.79), (4.87) and (4.89) to rewrite Sn in the form 

sn(.) ^ [ e»/m = If^/i.) - ^ I' t*'»'Im. (4.90) 

Taking into account (4.84) and (4.90), we have 



which implies 

d 

dx 



g{x) = {lCf,){x), IC := x-^l'^l-^ r^^N+i . ^^ g^^ 

^0 
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As /C is a bounded in i/(0, e) (1 < p < oo) Hardy operator, formulae (4.89) 
and (4.91) imply 

^{x)e 1/(0,8). (4.92) 

Prom (4.85), (4.86), and (4.92) follow relations (4.72). ■ 

P r o o f of Theorem 4.9. Using (4.56) and (4.57) one can see that (4.22) holds 
at Xq = 0. As go is bounded, the condition of Remark 4.2 is satisfied and 
identity (4.22) at = implies (4.23). Hence, the conditions of Proposition 
4.7 are fulfilled. Therefore, the GBDT generated by the parameter matrices 
A, S{0), and H(0) is a superposition of two GBDTs, where the first GBDT 
is generated by Ai, »Si(0), and ^i(O) (in addition to Proposition 4.7 see 
also (4.41)-(4.43) and Remark 4.6). As Si{0) > and qo is bounded in the 
neighborhood of zero, so the transformation of qo, which is generated by Ai, 
Si{0), \l/i(0), is bounded too. Denote this transformation by ql^\ 

The second GBDT in the superposition is determined by the x x x 
matrix A2 and by the matrix functions T22^{x), and 712(2;). As 5*21(0) = 
and Si{0) > formula (4.70) implies limj-^+o 7r2(a;) = ^'2(0). According to 
the second relation in (4.40) we have lim^^^+o ^22^ (a;) = 0, that is, we may 
put T2'^^(0) = 0. Thus, the second GBDT is generated by the parameter 
matrices ^2, ^22^(0) — 0> '^2(0) — ^2(0), which satisfy the identity 

^2^2-2^(0) - T2-/(0)^^ = ^2(0) J^2(0)*. 

If X = 1, then the potential ql^^ and the S'-node A2, T22^{0) = 0, and 
7r2(0) = ^1/2(0) satisfy the conditions of Lemma 4.10 and the statement of the 
theorem is true. 

If X > 1, then A2, ^22^(0) — ^2(0) satisfy the conditions of Lemma 
4.11, and the second GBDT is itself a superposition of GBDTs. Taking into 
account the block representation 

^W, 7r« e C (4.93) 

we derive that the transformation of gg^"* generated by A2, 0, and 7r2(0) is 
a superposition of the transformation generated by A^^\ 0, and 7r^^^(0), and 
of the transformation generated by A^2 \ 0) ^2^''(0)) respectively, where 
TT^^ is constructed analogously to the construction of tt2 in Proposition 4.7. 



Ao^ 



7^2 
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(The transformation of ^q^^ generated by A''^\ 0, and 7r(^)(0) is denoted by 
go^\) Moreover, for x > 2 the potential ql^^ and the S'-nodc ^2^\ 0, and 
712^'* (0) satisfy the conditions of Lemma 4.12, and according to Lemma 4.12 
this holds for further x— 3 steps too. At each step we obtain a new coefficient 
q^^ The absolute value of k'^^^ in the representation (4.71) of ^q'^'' (s > 1) 
equals s — 1, and the sign of k, changes at each step till we come to the final 
transformation. 

If X > 1, the final transformation is generated by 

4"^ = (A)..< e C, 0, 4-^0) = (0,0). 

Though we do not apply Lemma 4.12 at this last step directly, the proof of 
the representation of qo from Lemma 4.12 remains valid for qo — q^^^\ I 

5 Direct and inverse problems 

In this section we consider transformed systems constructed in Section 2, 
assuming that the initial systems are trivial, that is, either t> = or ,^ = 0. 
We use the notion of a minimal relization and some other notions from system 
theory, which are defined in Appendix A. 

5.1 Dirac systems 

5.1.1 Self-adjoint Dirac system 

First, consider the self-adjoint Dirac system constructed in Proposition 2.1 

^u{x,X) ^i{Xj + jV{x))u{x,X) (0<x<oo), 
We assume that 

5(0) = 4 > 0, A-A*^ m(o)jn(o)*, v{x) = o, (5.2) 

where v is the potential of the initial Dirac system (2.1). Partition 11 into 
two n X p blocks. Using (2.20), where y = 0, we derive 

n(x) = [*i(x) ^2{x)] = [e-^^^$i(0) e^^^*2(0)]. (5.3) 



V 
V* 



■ (5.1) 
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It follows from (2.7), (2.26), and (5.2) that 

v{x) = -2i$i(0)*e"^*5(a;)-^e^^^$2(0), 

S{x) ^In + [ U{t)U{t)*dt > 0. 
Jo 



(5.4) 
(5.5) 



As >S'(0) = /„ is fixed, parameter matrices in (5.3)-(5.5) are A and n(0) or, 
equivalently. A, $i(0), and $2(0). 

Definition 5.1 [27] The potentials v of the form (5.4) are called pseudo- 
exponential and the class of such potentials is denoted by PE. 

If V is locally summable on [0, cxd), there is a unique Weyl function of system 
(5.1) (see, for instance, [61]), which is defined in the following way. 

Definition 5.2 A holomorphic p x p matrix function ip such that 



[ Ip iifiXy ] Ku{x,Xyu{x,X)K* 



-i(fi{X) 



AeC+, u{Q,X) = Im, 



Ip Ip 
Ip Ip 



dx < 00, (5.6) 
(5.7) 



is called a Weyl function of Dirac-type system (5.1) on [0, 00). 

Moreover, the function (/?, which satisfies (5.6), is unique even without the 
analyticity requirement. 

Now, we shall again consider the case v G PE. As = 0, we have 
^(a;. A) = exp(ia;Aj) in (2.18), and so formula (2.18) takes the form 



u{x, X) = wa{x, A)e^^^^WA(0, A)-^ 



(5.8) 



By Proposition 6.2 [32] for every A e C, excluding a finite set of points, there 
is a limit 

\im wa{x,X) = fA^X). (5.9) 



If the equality 



K* 



-i(p{X) 



wa{0,X) 



j-p 




c(A) 



(5.10) 
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holds for somepxp matrices c(A) and <^(A) (A e C+), then formulas (5.8) and 
(5.9), and exponential decay of e*^^ imply inequality (5.6). As K* = K~^, it 
is easy to see that (5.10) is everywhere, excluding a finite set, equivalent to 
the relation 



Taking into account definitions (2.11) and (5.7) of wa and K, respectively, 
and formula (A. 4), one gets the following result (see [27] and Theorem 5.1 



Theorem 5.3 Let v G PE. Then the Weyl function of the Dirac system 
(5.1) admits realization 



ip{X) = Hp + 2$2(0)*(A7„ - i)-'$i(0), A:^A- i$i(0)($i(0) + $2(0))*. 



Note that the Weyl function 99 is a Herglotz function and admits represen- 
tation 



(^(A)=i. + |_" (^-^)dr(i), (u^u*,r^). (5.13) 



Here r is the spectral function of system (5.1). See Theorem 4.3 [27], where 
this fact is proved and an explicit expression for r in terms of the parameter 
matrices is given. The inverse problem to recover v G PE from r is solved 
explicitly in Theorem 4.5 [27]. (The corresponding bound states are con- 
structed in [31].) The inverse problem to recover v from the spectral density 
for the particular case a (A) C C+ was treated earlier in [5]. 

The left reflection coefficient for system (5.1), where v G PE, is expressed 
via the Weyl function by the formula (see p. 33 in [32]): 




-1 



(5.11) 



in [32]). 



(5.12) 



RL{X) = -{lp + t^{X)){lp-z^{X)) 



1 



(5.14) 



Recall that ip has properties 



%(A) > (AG C+), lim (^(A) = Hp. 



(5.15) 



46 



It follows from (5.14) and (5.15) that 

\\Rl{X)\\ < 1 (A e C+), lim Rl{\) = 0. (5.16) 

A— >oo 

To recover system (5.1) from its Weyl function we should recover param- 
eter matrices A and $fc(0) {k = 1, 2). First, we recover Rl via formula (5.14). 
By the second relation in (5.16) Rl is a strictly proper matrix function and 
admits a minimal realization 

RL(X)^C(XIn-A)-'B, (5.17) 

where C is a p x n matrix and B is an n xp matrix. According to the first rela- 
tion in (5.16), Rl admiting the minimal realization (5.17) is also contractive. 
Hence, by Theorems 21.1.1 and 21.1.3 in [LR] the Riccati equation 

XC*CX -i{AX - XA*) + BB* =0 (5.18) 

has a positive solution X > 0. Put 

X-^AX^ +iX-^BB*X-^, $i(0) = A'-^B, $2(0) = -iA'^C*. 

(5.19) 

The second relation in (5.2) is immediate from (5.18) and (5.19), that 

is, A and $^(0) {k = 1,2) are some parameter matrices and we can apply 
Theorem 5.3. It easily follows (see Section 9 in [32]) that (/? given by (5.12) 
and (5.19) coincides with the ip from which Rl was recovered. 

Theorem 5.4 [32] Let ip he a rationalpxp matrix function, which satisfies 
(5.15). Then p is the Weyl function of some system (5.1), where v G PE. 
To recovery we take a minimal realization (5.17) of Rl given by (5.14) and 
define parameter matrices by (5.19), where X > Q is a solution of the Riccati 
equation (5.18). After that we apply formulas (5.3)-(5.5). 

By Theorem 5.4 [32] all the potentials v G PE admit representation (5.4), 
where parameter matrices have additional properties 

n— 1 n— 1 

a{A) C C+, span |J Im yl'^$i(0) = C", span |J Im ^^=$2(0) = C", 

A;=0 fc=0 

(5.20) 

and Im means image. Hence, without loss of generality we can assume (5.20). 
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5.1.2 Potentials with singularities 

Now, consider a more general case of Dirac systems, including systems with 
singularities. Namely, let 

Ip, 



-Ip, 



s{o) = So, ASo - SoA* = m(o)jn*, (5.21) 

where pi, P2 > 0, pi+p2 — rn. Formula (5.5) takes the form 

S{x) = So+ f U{t)U{t)*dt > 0. (5.22) 
Jo 

As before, 11, wa, and v are defined via formulas (5.3), (2.11), and (5.4), 
respectively, and we assume that (5.20) holds. Here $i, $2, and v are n xpi, 
n X p2, and pi x p2 matrix functions, respectively. 

Definition 5.5 [22] The potentials v of the form (5.4), where (5.20) and 
(5.21) hold, and S is given by (5.22), are called generalized pseudo-exponential 
and the class of such potentials is denoted by GPE. 

Dirac system with v G GPE was treated on the whole axis in [62], and 
results on supertransparent potentials and soliton-positon interactions were 
obtained. In particular, it was shown in the proof of Theorem 2.1 from [62] 
that for some G M 

S{x) > for x> xq. (5.23) 

If 5*0 ^ 0, then det S{xk) = in some point (or points) on ]R_|_ U {0}, 
which means that the GPE class includes potentials with singularities. Nev- 
ertheless, as S{x) admits continuation meromorphic in x, the determinant 
det S{x) turns into zero only in a finite number of points on U {0}. The 
proof that the matrix function 

u{x,X) =WA{x,X)e''-'^^ (5.24) 

satisfies (5.1) and is nondegenerate (excluding a finite number of values of 
A and zeros of S{x)) remains true for Dirac system with j given by (5.21). 
Moreover, u as well as S is meromorphic in x. Therefore, we call u a fun- 
damental solution of (5.1), which agrees with the standard [45] requirement 
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for the fundamental and special solutions in the case of singularities to be 
defined by the same formula on the whole domain. 

The asymptotics of wa in (5.24) under condition (5.20) is given in the 
next proposition (see Proposition 3.1 in [22] or Theorem 3.1 in [62]). 

Proposition 5.6 Let A, S{0), and $^(0) {k = 1,2) satisfy (5.20) and 
(5.21). Then there is a limit 



and we have 
lim w^(x, A) — 



uj — lim ( e' 

x—*oo \ 







■ixA 



S{x)e 



ixA 



r 



>o, 



(5.25) 



X(A) 



X(A) := Ip, + i<^2{0)MA - Xlnr'^O). 

(5.26) 



Transmission and reflection coefficients for system (5.1) are defined [22,32,81] 
in terms of special solutions of (5.1), that is, m x pi and m x p2 solutions y 
and Z, respectively, which we determine by the relations 



y{x,X)=e 
Put 



ixX 







+ o(l) (x^oo), Z{0,X) 







'P2 



A e 



Tl(x) := yiio, x)-\ Rl(x) := :^^2(o. A) y,{o, A) 



-1 



Rr{X) :-- 
Tr{X) := 



( lim 

\x—*oo 



e-^^^Zi(x,A))( lim e*^^Z2(x 



^,A))"\ 



ixX ■ 



f lim e'^^Z2ix 

Vx— >oo 



;,A))"\ 



(5.27) 

(5.28) 
(5.29) 

(5.30) 



Here the pi x pi matrix yi and the pi x p2 matrix are upper and lower 
blocks of y, respectively. Analogously, the pi x p2 matrix Zi and the p2 x p2 
matrix Z2 are upper and lower blocks of Z, respectively. The functions Tl 
and Tr are called the left and the right transmission coefficients and Rl and 
Rr are called the left and the right reflection coefficients, respectively. 
Using (5.24) and (5.26) we get the result 
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Theorem 5.7 [22] Let v e GPE. Then the transmission and reflection 
coefficients are given by the formulas 





= Jp, +^<l>i(0)*5o- 




-A/„)-'$i(0), 


(5.31) 


Rl{X) 






r^$i(o), 


(5.32) 


Tb{X) 


= /p,+Z$2(0)*5o 


■\e- 


-\In)-\ln-SoCu)^2{0), 


(5.33) 




= ^<l>^{OrS,\e- 




1-1(4 -5o^)$2(o) 












(5.34) 



where 9 -.^ A - i^i{0)^i{0)*S^^ . 

The inverse problem to recover v from Rl is solved in Theorem 4.1 [22]. 

Theorem 5.8 Let TZ be a strictly proper rational p2 x pi matrix function. 
Then TZ is the left reflection coefficient of a system {5.1) with v e GPE if 
and only if TZ is contractive on M. If TZ satisfies this condition, then v can 
be uniquely recovered from TZ in two steps, that is, steps (i) and (ii) below. 

(i) First, take a minimal realization 

TZ{\) = C{\In - A)-^B. (5.35) 
Then there is a unique solution X of the Riccati equation 

i{XA- A*X)^C*C + XBB*X (5.36) 

such that 

a{A + iBE*X) ^C_, X = X\ (5.37) 

Moreover, det X ^Q. 

(ii) Next, recover parameter matrices by the equalities 

A^A + iBB*X, So = X-\ ^i{0) = B, $2(0) = -i^oC*, (5.38) 
and recovery via {5.3), {5.4), and {5.22). 
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5.1.3 Skew-self-adjoint Dirac system 



In this subsubsection we consider system 
-^u{x, A) = {iXj + jV{x))u{x, A) (0 < X < oo), °^ 



, (5.39) 



v{x) = 2$i(0)*e^"^*5(x)-^e^"^$2(0), S{x) ^In+ f U{t)3Tl{t)*dt, 

Jo 

(5.40) 

where j has the form (1.7), 11 is given by (5.3), and the identity A — A* — 
in(0)n(0)* is fulfilled. That is, we consider GBDT transformation of the 
trivial initial system, where v = 0. According to (2.23) and to the second 
relation in (5.40) the inequality S{x) > holds (see inequalities (1.6) and 
(1.7) in [28]). Thus, the right-hand side of the first equality in (5.40) is well- 
defined. The class of potentials v of the form (5.40) is denoted by PE2. By 
Proposition 1.4 in [28] for each v G PE2 there is some M > such that 

sup \\v{x)\\ < M. (5.41) 

xE[0, 00) 

It follows from Proposition 2.2 that the fundamental solution u{x, A) of 
system (5.39) (normalized by u{0, A) = Im) is given by formula (5.8), where 
wa has the form (2.24). 

Analogously to Definition 5.2 one can define Weyl functions of the skew- 
self-adjoint Dirac system [16,28,52,68]. 

Definition 5.9 Let system {5.39), where \\v{x)^ is bounded on all the finite 
intervals [0, I], be given. Then, a holomorphic p x p matrix function (p such 
that 



POD 

/ [ifiXy Ip]u{x,Xruix,\) 
Jo 



V.(A) 



dx < 00, 9A < -Ml (5.42) 



is called a Weyl function of {5.39). 

There is a unique Weyl function of system (5.39) with a bounded on [0, 00) 
potential [52,68]. If (5.41) holds, we can put in (5.42) Mi = M. For a 
generalization of the Definition 5.9 see [16,23,53]. 

Direct and inverse problems for v G PE2 are solved explicitly in Theorems 
2.1 and 2.3 from [28]. It is done in a similar to the self-adjoint case way. 
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Theorem 5.10 [28] Let v e PE2. Then system {5.39) has a unique Weyl 

function ip, which satisfies {5.42) on the lower semi-plane C_, a finite number 
of poles excluded, and this function is given by the formula 

(^(A) =i$i(0)*(A/„-i)-^$2(0), i:= A-i$2(0)$2(0)*. (5.43) 

Theorem 5.11 [28] Let ip he a strictly proper rational p^p matrix function. 
Then ip is the Weyl function of the skew- self- adjoint Dirac system {5.4-2), 
where v e PE2. To recover v take a minimal realization 

ip{X) = C{XIn - A)-'B. (5.44) 

There is a positive solution X > of the Riccati equation 

XCXX + i{AX - XA*) - BB* = 0. (5.45) 

After putting 

$1(0) = iX'^C\ $2(0) = X-^B, A = X-^AX^+iX-^BB*X-^, (5.46) 
the potential v is recovered by formula {5.4O). 

5.2 AT- wave equation 

Here we consider GBDT transformation of the auxihary system for the N- 
wave equation 

^u{x,X)^{iXD-[D,C{x)])u{x,X) {0<x<oo), f = f ■ (5-47) 
We consider the subcase of (2.27), where B = Im, and assume that 

D = diag {di, ^2, ■ ■ ■ , dm}, di > d2 > . . . > dm > 0- (5.48) 

Starting from the initial system (2.27) with ^ = and taking into account 
B — Im, rewrite (2.31) as: 

^= n*5-^n, (5.49) 

where 

AS{o) - s{o)A* = in(o)n(o)*, = -iAUD, = n^n*. (5.50) 
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We require 

S(0,A) = 7^, S{0)>0. (5.51) 

As we have u{x, A) = exp{ixXD) for the fundamental solution of the initial 
system, formula (2.18) takes the form 

u{x, A) = wa{x, A)e*^^^WA(0, X)-\ (5.52) 

where wa has the form (2.24) and AS - SA* = inn*. It follows from (2.24) 
and AS - SA* = iUU* that 

wa{x, XywA{x, X)^Im- i{X - X)U{xy{A* - XIn)-'S{x)-\A - A7„)-^n(x). 

(5.53) 

By the last relations in (5.50) and (5.51) we get S{x) > for x > 0, that 
is, wa{x,X) is well-defined for all a; > and all A ^ (^{A). If A e C_, the 
inequality S{x) > and formula (5.53) imply 

WAix, X)*wa{x, A) < Im, (5.54) 
i{X - X)U{x)*{A* - XI^y^Sixy^A - A7„)-^n(x) < Im- (5.55) 

Proposition 5.12 Let parameter matrices A, S{0), andll{0) satisfy condi- 
tions AS{0) - S{0)A* = in(0)n(0)* and S{0) > 0. Then ^ determined by 
{5.49) and {5.50) is hounded on [0, 00). 

Proof. It follows from (5.50) that li has the form 

n(x) = exp ( - ... eKY>{- idmxA)fm , (5.56) 

fkeC^ {l<k<m). 

One can easily see that span^^^lA — XIn)~^fk ^ fk foi' any open domain O. 
Hence, by (5.55) and (5.56) we have 

sup ||^(x)-^e-^"''=^^/fe|| < 00 (1 < A; < m). (5.57) 

a;GM+ 

Now, the boundedness of ^ given by (5.49) is immediate. ■ 

The class of potentials ^ of the form (5.49), which are generated by 
parameter matrices satisfying conditions of Proposition 5.12, is denoted by 
PE3. 
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Definition 5.13 [59, 73] A Weyl function of system (5.47) is an m x m 

matrix function 99, such that for some M > it is analytic in a lower semi- 
plane SJA < —M, and the inequalities 

sup \\u{x, X)ip{X) exp {—ixXD)\\ < 00 (5.58) 

x<l, 3iA<-M 

hold for all I < 00. 

In our case ^ is bounded, and so (see [59]) the system (5.47) has a unique 
Weyl function ip normahzed by the condition 

^,,(A) = 1 for /c = ^,,(A) = for > s. (5.59) 

This Weyl function satisfies for some £ > the inequahty 

00 

/ (exp (i.-XD))mM., A)-5(x, A)^(A) axp (x(-.AB - elj) d. < oo. 



(5.60) 

Here, we do not require (5.59). Our next theorem is immediate from (5.52) 
and (5.54). 

Theorem 5.14 A Weyl function of system {5.4'T}, where ^ e PE3, is given 
by the formula 

ip{\) = w;^(0, A) = 7^ - in(0)*^(0)-^(^ - A/„)-^n(0). (5.61) 

Notice that in view of (5.53) and (5.54) the matrix function (^(A) given by 
(5.61) has following properties: 

(^(A)(^(A)* = 7^, (^(A)V(A)<7^ (AeC_), lim <^(A) = 7^. (5.62) 

A— »oo 

In particular, this matrix function admits a minimal realization 

ip{\) + C{\In - A)-^B. (5.63) 

It follows from the second relation in (5.62) that (y{A) G C_|_. Hence, there 
is a unique and positive solution ySo > of the identity 

ASo - SoA* = iBE\ (5.64) 
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Matrix functions admitting realization (5.62) satisfy conditions of The- 
orem 1.2 from [59] and there is at most one solution of the corresponding 
inverse problem. This solution is given in Theorem 5.6 from [73]. 

Theorem 5.15 [73] Let apxp rational matrix function ip satisfy conditions 
{5.62). Then (f is a Weyl function of the unique system (5.47), where ^ e 
PE3. To recover ^ take a minimal realization {5.63) of (p, recover Sq from 
{5.64), and put 

A^A, S{0) = So, n(0) = B. (5.65) 

Then^ is generated by the parameter matrices A, S{0), andll{0) via formulas 
{5.49) and {5.50). 

Finally, let parameter matrices A, 5(0, 0) > 0, and 11(0, 0) satisfy (2.36), 
where B — 7^, and generate via (5.49) the solution ^ of the N-w&ye equation 
(2.32). The corresponding evolution of the Weyl function is given by the 
formula 

^{t, A) = wa{0, t, X) = Im- ^^(0, t)*5(0, t)-\A - A4)-^n(0, t). (5.66) 

Here Ht — —iAUD, St — TIDW , and formula (5.66) holds on the interval 
t e [0,e), where S{0,t) > 0. 
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A Mathematical system theory 

We present here some basic results and notions from mathematical system 
theory of rational matrix functions that are used in our review. This material 
has its roots in Kalman theory [37], and can be found in various books (see, 
for instance, [8,17]). See also interesting historical remarks in [39]. 

The rational matrix functions appearing in the article are proper, that is, 
analytic at infinity. Such an m2 x mi matrix function W can be represented 
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in the form 

W(X)^V + C(XIn-A)-^B, (A.l) 

where ^ is a square matrix of some order n, the matrices B and C are of sizes 
nxmi and 777-2 x n, respectively, and V — l^(oo). The representation (A.l) 
is called a realization of W, and the number ord(^) (order of the matrix A) 
is called the state space dimension of the realization. 

The realization (A.l) is said to be minimal if its state space dimension n 
is minimal among all possible realizations of W. This minimal n is called the 
McMillan degree of W. The realization (A.l) of W is minimal if and only if 

n— 1 n— 1 

span y Im ^'^B = C^, span |J Im (X)*^C* = C", n = ord(^), (A.2) 

k=0 k=0 

where Im is image. If for a pair of matrices A, B the first equality in (A.2) 
holds, then the pair A, B is called controllable or a full range. If the second 
equality in (A.2) is fulfilled, then C, A is said to be observable. If a pair A, 
B is full range, and /C is an mi x n matrix, then the pair A — BfC, B is also 
full range. 

Minimal realizations are unique up to a basis transformation, that is, if 
(A.l) is a minimal realization of W and if W{X) = T) + C{XIn — A)~^B is 
a second minimal realization of W, then there exists an invertible matrix S 
such that 

A--SAS-\ B^SB, C^CS-\ (A.3) 

In this case, (A.3) is called a similarity transformation. 

Finally, if is a square matrix and V — (mi = m2), then 
admits representation 

W{X)-' = Im^-CiXIn-A'^yB, A'^=A-BC. (A.4) 
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